Critical Current Enhancement due to an Electric Field in a Granular d-Wave 

Superconductor 
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We study the effects of an electric-field in the transport properties of bulk granular superconductors 
with different kinds of disorder. We find that for a d-wave granular superconductor with random 
7r-junctions the critical current always increases after applying a strong electric field, regardless of 
the polarity of the field. This result plus a change in the voltage as a function of the electric field 
are in good agreement with experimental results in ceramic high T c superconductors. 



PACS numbers: 74.50. +r, 74.60. Jg, 74.80. Bj 

Recent experiments in ceramic high-T c superconduc- 
tors (HTCS) have found an enhancement in the critical 
current when applying an electric field E through an in- 
sulating layer jl]-[|] . Previous studies of the electric field 
effects (EFE) in superconducting films have attributed 
the changes in the critical current to variations in the 
charge density or to a redistribution of carriers, which ap- 
pear at the surface layer with depths of the order of the 
electrostatic screening length <1e (in HTCS, rig « 5 A) 
0. Several experiments in ultrathin YF^CusCv-z films 
(5 — 10 nm thick) have found that E can affect T c and 
the IV characteristics in good agreement with this 
picture ||. In this case, there is either an enhancement 
or a depletion of the critical current depending on the 
polarity of E 0] . The surprising observation in [|l]-[| of a 
strong EFE in bulk ceramic HTCS (1.5 mm thick), how- 
ever, can not be explained by a surface effect. Moreover, 
for high enough electric fields, the critical current always 
increases regardless of the polarity of the field . 

Rakhmanov and Rozhkov have shown that an elec- 
tric field can induce a change in the critical currents 
of the Josephson junctions present in granular samples. 
However, in their model the critical current either in- 
creases or decreases depending on the sign of E. Re- 
cently, Sergeenkov and Jose Q have proposed that an 
electric field applied to a granular superconductor can 
produce a magneto-electric like effect, which could be in- 
directly related to the behavior of the critical current in 
[Q-Q, but no comparison with the experimental results 
was given. 

Granular superconductors are usually described as a 
random network of superconducting grains coupled by 
Josephson weak links [pl[L0|. In the HTCS ceramics 
several experimental groups have found a paramagnetic 
Meissner effect (PME) at low mag netic fields [JO). Sigrist 
and Rice ]T2] proposed that this effect could be a conse- 
quence of the intrinsic unconventional pairing symmetry 



of the HTCS of d x 2_ y 2 type p3| . Depending on the rela- 
tive orientation of the superconducting grains, it is possi- 
ble to have weak links with negative Josephson coupling 
(^-junctions )fl2|,[l3| which, according to ]lj]Ji^| , give rise 
to the PME |14[|. In this paper we will show that the pres- 
ence of 7r-junctions in ceramic samples also explains the 
unusual electric field effects observed in |l] [|] . 

We consider a 3-D cubic network of superconducting 
grains |],[l(| at n = (n x ,n y ,n z ) and with unit vectors 
pb = x,y, z. The current //i(n) between two grains n, and 
n + /i is given by the sum of the Josephson supercurrent 
plus a dissipative Ohmic current: 
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Here # M (n) = #(n + fi) — 6(n) — /^(n, t) is the gauge 
invariant phase difference, with 0{n) the superconduct- 
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ing phase in each grain, and ^4 M (n, t) = |^ J" 
($0 = h/2e). The critical current of each junction is 7° „ 
and R is the normal state tunneling resistance between 
grains. Together with the conditions of current conser- 
vation, J2/j,[^i^( n ) ~ ^( n — A)] = 0' tn i s determines the 
dynamical equations for the Josephson network. We con- 
sider periodic boundary conditions (PBC) in a network 
with N x N x N grains. 

When an electric field E is applied in the z direction, 
the z component of the vector potential is given by: 



2ncd „ 

A z {n, t) = A z (n, 0) — — Et, 



(2) 



with d the intergrain distance or junction thickness. 
This results in a high-frequency alternating supercur- 
rent in the z-direction due to the ac Josephson effect 
p5[ . In addition, we consider that the sample is driven 
by an external current density I ex t along the y direc- 
tion. Therefore, the vector potential term is A M (n, t) = 
—Sfj,, z uJEt — Sfj^yCtylt), where the electric field frequency 
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is uje — 2ircEd/$o, and we consider that no external 
magnetic field is applied. (The experiments had zero 
magnetic field The external current density with 

PBC determines the dynamics of a y (t) as 
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The average voltage per junction induced by the driving 
current is then obtained by V = 

Furthermore, we consider that the applied electric field 
is screened inside the grains and acting only in the insu- 
lating intergrain region of the junctions, of typical thick- 
ness d ~ 10 — 20A. We also neglect the effects of inter- 
grain capacitance Cj and intragrain capacitance C g . In 
general, the effect of the capacitances is to screen the ap- 
plied field E ext inside the sample within an overall screen- 
ing length \e ~ (Ci/Cg) 1 / 2 . Since C g <C Cj, \e is very 
large and we can consider that the internal field acting 
in the intergrain junctions is = E w a p E ext with 
the polarizability a p < 1. We also neglect the possible 
E dependence of the critical currents of the junctions. 
As shown in JtJ this assumption may result in an in- 
crease or decrease of J° depending on the sign of E in . 
We also neglect screening current effects (i.e. finite self- 
inductances). The self-induced magnetic fields are very 
important for the description of the critical state [l6| at 
large magnetic fields, and for the PME at low magnetic 
fields 0. The experiments of 0-§ were done at zero 
magnetic field and finite electric fields. Thus the most 
important approximation we make is to neglect the ca- 
pacitances (i.e. the screening of E). After making all 
these physical approximations, we will now focus on the 
collective aspects of the ac Josephson effect induced by 
the electric field. In this model, the electric field scale is 
given by E a = RI a (T = 0)/d= 7rA s (0)/2ed, with A s (0) 
the superconducting energy gap; for the YBaCuO ceram- 
ics we have A s (0) « 20meV, which gives Eq ~ 3QMV/m, 
i.e. in the same range as the fields used in Jl|-Q]- 

We assume that the effect of disorder in a granular 
superconductor at zero magnetic field mainly modifies 
the magnitude of the critical currents. In s-wave super- 
conductors the sign of the Josephson coupling is always 
positive. In d-wave ceramics the sign of I® „ is expected 
to vary randomly depending on the relative spatial orien- 
tation of the grains. Therefore, we consider two models 
of disorder: (i) Granular s-wave superconductor (GsS): 
We consider 7° a random variable with a uniform dis- 
tribution in the interval [Io(l — A c ),io(l + A c )] with 
(J„ jA1 ) = Io and A c < 1. (ii) Granular d-wave super- 
conductors (GdS): We consider that there is a random 
concentration c of 7r-junctions with 7^ „ = — Io, and 
1° = Io with concentration 1 — c 17-19). The latter 
model was used to explain the history-dependent param- 
agnetic effect the anomalous microwave absorption 
Jl8[ and the non-linear susceptibility and glassy behavior 



p9| observed experimentally in ceramic HTCS pij ]. 

We have performed numerical simulations of systems 
of sizes N — 8, 16 with a numerical integration of the 
dynamical equations with time step At = O.Itj, with 
tj = ^o/2ttcRIq, and we took 5 x 10 4 integration steps. 
We calculated the current-voltage characteristics (IV) for 
different amounts of disorder and electric fields, averaging 
over 20 realizations of disorder in each case. 

In the absence of disorder, the IV curves are unaffected 
by the electric field. In a perfect cubic network, a fi- 
nite electric field induces an ac supercurrent Iq sm(u> Et) 
along the z axis only, and therefore the IV curves in the 
xy plane are not affected by the value of E. When the 
disorder is small, the amplitude of the ac supercurrents 
in the z direction is random and, due to current conser- 
vation in each node of the network, this induces small 
ac-currents in the xy planes with random amplitudes 
~ A c Jo sin(a>.Et). Adding a small ac current to a Joseph- 
son junction reduces its effective dc critical current [ [l5|| . 
In fact, this is the effect we observe in the GsS model. 
In Fig. 1(a) we show the IV curve for A c = 0.6. There 
we see that after applying an electric field E the whole 
IV curve shifts to lower values of the current, the appar- 
ent critical current decreases, and therefore the voltage 
change AV = V(E) — V(0) is positive for a given current 
I > l c . However, this is in the opposite direction of the 
experimental results of El-0], where an increase in the 
critical current was observed. 
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FIG. 1. Current- voltage characteristics before and after ap- 
plying an electric field E. (a) For a granular s-wave supercon- 
ductor (GsS), with A c = 0.6. (b) For a granular d- wave super- 
conductor (GdS) with concentration c = 0.5 of 7r-junctions. 
Voltages are normalized by NRIq and currents by N 2 Io and 
lattice sizes are 16 x 16 x 16. 

Let us then consider the case of a granular d-wave 
superconductor. The presence of randomly distributed 
negative and positive critical currents changes the previ- 
ous simplistic picture, since now frustration is introduced 
in each loop containing an odd number of 7r-junctions 
@,|3Jl9|. In Fig. 1(b) we show the IV curve for the GdS 
model with a concentration c = 0.5 of 7r-junctions. We 
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note that after applying the electric field the IV curve 
is shifted towards higher currents, i.e. the "apparent" 
critical current increases. If we change E — > —E the IV 
curves overlap showing that the effect is independent of 
the polarity of the electric field. A similar effect was seen 
in the experiments of ||: the IV curves shift upwards 
for large E and almost overlap when changing the po- 
larity of the field [^0| . This result is surprising since in 
the usual electric field effect observed in films a change 
in the polarity of the field drastically changes the sign 
of the shift of the IV curves 0. We also see in Fig. 1(b) 
that the EFE is stronger near the critical current, and 
for large currents the change in voltage AV is smaller, 
which was also seen in experiment Jl]-^| . If the amount of 
disorder in the GsS is such that A c > 1, then a fraction 
(A c — 1)/2A C of the critical currents will be negative. 
This case is unrealistic for s-wave superconductors, but 
corresponds to a d-wave granular superconductor with a 
small concentration of 7r-j unctions. In fact, we also find 
an increase of the critical current for A r > 1. 




FIG. 2. Changes in voltages as a function of electric field 
for different kinds of disorder. AV/V = (V(E)~V(0))/V(0). 
The results for the GsS model are for I = 0.97o and for the 
GdS model, I = 0.5/o- Electric fields are normalized by -Co- 
in Fig. 2 we study the relative change in voltage 
AV/V(0) for a fixed bias current (above the critical cur- 
rent) as a function of E. We see that for the GsS model 
AV is positive, and only when negative critical currents 
are allowed (A c > 1), there is a decrease in the voltage, 
i.e. AV < 0. The effect is stronger in the GdS model 
with c = 0.5 concentration of 7r-junctions. An effect of 
the same order, i.e. near a 50% decrease in voltage, and 
in the same scale of electric fields was observed in experi- 
ment. In Fig. 3 we show the dependence of the EFE with 
different amounts of disorder. For a given bias current 
we calculate AV/Vq after applying a large electric field 
of E = 0.8Eq. We see in Fig. 3(a) the voltage change in 
the GsS model as a function of A c . It is clear that only 
when A c > 1 the AV can be negative. In Fig. 3(b) we 
study the same situation for the GdS model as function 



of the concentration of the 7r-junctions. We obtain that 
a small amount of 7r-junctions is enough to show a strong 
decrease in voltage at large electric fields. We see that 
there is a very rapid change in AV/Vq close to c = 0.1. 
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FIG. 3. Effects of disorder in the voltage change after ap- 
plying an electric field E. (a) GsS model, AV/Vq vs. A c , 
I = 0.9J . (b) GdS model, AV/V vs. c, I = 0.5I . All the 
results were obtained with a field E = 0.8i?o. 

It is interesting to note that there is another type of dis- 
order that can give a negative AV/Vq. In the presence of 
external magnetic fields, the vector potential has a com- 
ponent A^(n). For strong magnetic fields, Ajf(n) can 
be taken as a random variable in the interval [0, 2% Ah], 
with the case A H = 1 corresponding to the gauge glass 
model J2l|] . We have found that only when Ah > 0.5 we 
can obtain a decrease in the relative voltage, AV/Vq < 0, 
similar to the results obtained with the GdS model. The 
case Ah > 0.5 is when the effective couplings can take 
negative values. This case, even when interesting in it- 
self, is unrealistic for this problem since it corresponds to 
ceramics at very large magnetic fields, while the experi- 
ments had zero magnetic field. 

There are other interesting aspects of the experimen- 
tal results. It was found that for low electric fields the 
critical current decreases and only for large fields it in- 
creases. A typical measurement shows that, when apply- 
ing a current near the critical current, AV/Vo increases 
as a function of E and, after reaching a maximum, it 
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falls and then becomes negative for large enough fields. 
Our random d-wave model for a granular superconductor 
also reproduces this effect in good agreement with exper- 
iment as we show in Fig. 4. We note that the value of 
AV/Vo at the maximum decreases when increasing the 
bias current, and only for large enough currents, the volt- 
age change AV/Vo is always negative and decreasing with 
E as previously shown in Fig. 2. The same dependence 
with the bias current has been observed in experiment 
(see for example Figs. 3 and Fig. 4 of Q| and Fig. 3 of 0). 
These results can be understood by looking at the shape 
of the IV curves at low electric fields. As we show in 
the inset of Fig. 4, for low E there is a crossing of the IV 
curves, which explains the behavior observed in AV/Vq. 
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FIG. 4. Change in voltage as a function of E for differ- 
ent driving currents for the GdS model with c = 0.5. Inset: 
Crossing of IV curves for different electric fields. 

The experiments on the electric field effects in ceramic 
HTSC were surprising for two reasons: (i) a large electric 
field effect in a bulk ceramic sample was unexpected [jjj, 
(ii) an increase in the apparent critical current as a func- 
tion of E was found to be independent of the polarity 
of E ||. Here we have qualitatively explained all these 
experimental features with a simplified model of a granu- 
lar superconductor, which has two basic ingredients: the 
ac Josephson effect induced by the electric field and the 
collective frustration effects due to 7r-junctions. There 
are many interesting open questions such as: the coexis- 
tence of this effect with the paramagnetic Meissner effect 
pl| , P^Jl7| , the time dependent glassy dynamics as a func- 
tion of electric and magnetic fields the effects in a 
chiral glass state [[l9| , and the effect of finite inductances 
and capacitances in a more realistic model. We expect 
that the results presented here will motivate further ex- 
perimental and theoretical studies of this problem. 
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